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O Abstract 

Several general results for the spectral determinant of the Schrodinger operator on metric 
Cn) graphs are reviewed. Then, a simple derivation for the (^-regularised spectral determinant is 

proposed, based on the Roth trace formula. Two types of boundary conditions are studied : 
^ functions continuous at the vertices and functions whose derivative is continuous at the ver- 

'""5 tices. The ^-regularised spectral determinant of the Schrodinger operator acting on functions 

with the most general boundary conditions is conjectured in conclusion. The relation to the 
Ihara, Bass and Bartholdi formulae obtained for combinatorial graphs is also discussed. 

r—n PACS numbers : 02.70.Hm ; 02. 10. Ox 
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^ 1 Introduction 

a 

^ Spectral determinants (i.e. (^-functions or L-functions) and trace formulae in graphs have 

I 1 attracted a lot of interest in the mathematical literature p2 l [M l [35 l [20l.[2Tl[3l lT5l[371 [2ll5l[T7 1 

^ [251 19] and in the physical literature as well [321 [Ml II [261 [IJ [JOl [III El [Ml [HI SSI [7] . In this 

^ article we will consider the case of the Laplace —A or the Schrodinger operator — A -|- V(x) on 

a metric graph. We denote by {E^} the spectrum of the operator —A -|- V{x), and introduce 
the spectral determinant by the formal definition 

(N 

O 
O 



X 



^(7) = det[7 - A + Vix)] = JJIt + En) , (1) 



where 7 is some "spectral parameter". This object was shown to be a very convenient quantity 
in order to study several questions in metric graphs, like various properties of the Brownian 
motion [lOl [HI [H [121 [7] or transport and magnetic properties of networks of metallic wires 
I32l[33l[ll[ll|. 

^ Because the operator —A -|- V{x) acts in a space of infinite dimension, the computation of 

the spectral determinant ([T]) requires in practice some regularisation. Starting from the trace of 
the Green function, i.e. the Laplace transform of the partition function Z(t) = X^n^"*^"' 



1 roc 
G(7) = V — = / dtZ 

^1 + En Jo 



it)e 



-7t 



furnishes a first possible regularisation, used in Refs. [32l [331 E [3 [13 [HI [3 HI]. Performing 
some integration with respect to the spectral parameter we obtain 

5GF(7) exp dV G(7') = S(7)/5(7o) (3) 
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For example, consider the Laplace operator on a line [0, L] with Neumann boundary conditions. 
The spectrum is in this case En = (nvr/L)^ with n G N and we obtain (7) = ^^sinh^\ - 
This procedure however leaves some arbitrary in the choice of the parameter 79. 

Another well known regularisation for determinants is the ("-regularisation whose starting 
point is to introduce the ("-function 

C(^, 7) = E + ^-)'' = / t^-' Z(t) , (4) 

i [s) Jo 

that we have expressed as a Mellin transform of the partition function, for convenience for the 
following discussion. Note that ((1,7) = ^(7). The ("-regularised determinant is then related 
to a derivative of the (-function : 

S^(7) = exp-^(0,7). (5) 

In this case no arbitrary is left within the calculation of the determinant. Coming back to 
the simple example considered above of the Laplace operator on a line [0, L] with Neumann 
boundary conditions, we get in this case 5^(7) = 2^/7sinhTy7L (derived below). 

An important remark is that various regularisations only differ by some 7-independent pref- 
actors. This explains why the choice of the regularisation has no consequence on the properties 
studied in Refs. j32l [33l [U [IHl [HI EJ El HSJ [7] since they are always related to derivatives of 
ln5'(7), with respect to the spectral parameter 7 or other parameters. In a recent work [44J 
some procedure was proposed in order to construct the spectral determinant of a graph by 
combinations of the determinants of subgraphs. In this case it is crucial to define precisely the 
prefactor of the spectral determinant. 

In a recent paper [T7], Friedlander has derived the (-regularised determinant for the Schro- 
dinger operator on a metric graph. One should also mention that the analysis of the regularised 
determinant of the Laplace operator for general boundary conditions, det'(— A) with the prime 
indicating exclusion of zero mode contribution is there is some, was the subject of the very 
recent work [19]. It is the purpose of the present article to propose another derivation of the 
(-regularised determinant for the Schrodinger operator det[7 — A + y(j;)]. Our approach is based 
on the Roth trace formula [34] and will allow simple extension of Friedlander's result to other 
choices of boundary conditions at the vertices. 

In the next section we set notations. In section [3] we mostly recall some known general results 
obtained by Desbois in Ref. [10] and needed for the following sections. The section |4] focuses 
on the case of functions continuous at the vertices and section [5] on the case of functions with 
derivative continuous at vertices. The (-regularisation of the spectral determinant is provided 
in section [H 



2 Metric graphs and Laplace operator 

A graph is a set of V vertices (here labelled with greek letters a, f3,...) linked by B bonds 
(denoted as (a/3),...). Each bond (a/3) is associated to two arcs (oriented bonds) a/3 and /3a (arc 
will be also labelled with roman indices a, b,....). For the arc a = af3, the reversed arc is denoted 
a = /3a. We introduce the adjacency matrix : Oq,^ = 1 if a is linked to /3 by a bond and 0^/3 = 
otherwise. We denote by rua = 'Ylji'^oiii the coordination number (valency) of the vertex. The 
graph is said to be a metric graph (or a quantum graph) if each bond is identified with a finite 
interval (0, lap) G M, where lap designates the length of the bond linking the vertices a and /3. 
In this case we may consider a scalar function ip{x), characterised by a set of B components 
ipapixap) with Xap S {OJa/s)- The Component is labelled by the arc in order to specify the 
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direction along which the coordinate is measured [this redundancy of the notation imphes the 
obvious relation ipapixap) = ippaixpa) = ipiSaiLis - Xais)]- 

Having introduced scalar functions living on the graph we may define the action of the 
Laplace operator A on such functions. Along a wire it acts as the usual second derivative 
(A^)q,^(x) = ^il^ai3{x)- At the vertices, the set of functions on which A acts must satisfy 
some boundary conditions in order to ensure self adjointness of the operator. Let us denote 
by ip{0) and tp'{0) the vectors of size 2B gathering the values taken by the components and 
the derivatives at the vertices V'a/3(0) and ^'0/3(0); respectively. The most general boundary 
conditions ensuring self adjointness of A are of the form 

CVW + DV^'W = 0, (6) 

where the 2B x 2B matrices C and D satisfy (i) CD^i = DC'< . (ii) The 2B x AB matrix {C, D) 
has maximal rank (equal to 2B) [24J. Note that, characterising what arc is connected to what 
other arc, these two matrices encode all the information on the topology of the graph. 



3 General results for the spectral determinant 



In this section we mostly recall some results obtained by Desbois [10] for the spectral 
determinant of the Schrodinger operator for general boundary conditions. We compute the 



spectral determinant by constructing the Green function Q{x,y) = {x\ 



y-D-i+V{x) 
2 • 



\y), where 



T)x = gj: — iA{x) is the covariant derivative. Note that the replacement of A by is motivated 
by physical considerations [311 [U [39] . It might also be useful in order to study winding proper- 
ties of Brownian curves in the graph [42^ [7] . It will affect calculations in a very simple manner 
through the introduction of additional (magnetic) phases. In the following it will be understood 
that the 2i?-vector ip'{0) of Eq. ^ gathers the covariant derivatives. 



a 







ai ^ 








Figure 1: The vertex a. (a) For general boundary conditions, matching of the Green function at 
the vertex is described through two matrices x rua denoted C^"-* and D^°'^ (the blocks of the 
matrices C and D). One must assign a value of the Green function Qa to each arc a. (b) For 
continuous boundary condition at the vertex, the matrices C^"^ and D^"^ take the form (29). In 
this case one can assign to each vertex a a single "vertex variable" Qa- 



3.1 Arc determinant (1) 

Let us consider two points x and y belonging to the two arcs a and b, respectively. On 
the arc a we use the set of independent solutions fa and fa of the differential equation [7 — 
^ + Va{x)]f{x) = such that /a(0) = 1 and fai}a) = (see appendix|A). The Green function 
depends on two coordinates x and y and therefore on must specify by two indices (here a and 
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5) to which arcs they belong to : 

Qa,b{Xa,yb) =Qa fa{Xa) 6^^°"" + Qa UXa) e'^^^^ 



Wa 



/a(max {xa,ya)) /a(max {Xa, Va)) , 



(7) 



where Qa and Qg, denote the values of the Green function at the two ends of the bond (a) (see 
figurejlja). We have introduce the notation (^(a),(b) = ^a,b + ^a,b- Note that the dependence of Qa 
and Qa in the coordinate y is implicit. Aa = A{x) for x S (a) is the (constant) vector potential 
on the bondj^ and 9a = Aala the magnetic flux along the wire. Wa = —f'a{la) is the Wronski 
determinant of the two linearly independent solutions, Eq. (69). All the values of the Green 



function at the vertices Qa^bi^^Vb) = Qa (figurejlja) are gathered in the vector ^(0) of size 2B 
and all covariant derivatives iJ^xQ)a,b{Q-,yb) in the vector ^'(0). We obtain 



Q'{{)) = MQ{{)) + 



\ 





h{yb)e-'^''y' 



V 



(8) 



where the matrix M couples an arc a to itself and to the reversed arc a : 

Ma,b = Sa,bf'a{0)-6rfL{la)e''-. 



(9) 



Noticing that f'a{la) = f'ai^a) G IK (see appendix [A|) and 9a = —9a shows that this matrix is 
Hermitian. Imposing boundary conditions at vertices, C Q{0) + D Q'{0) = 0, we obtain all values 
Qa, that could be reinjected into Eq. ([T]). The Green function at coinciding points reads : 



QaA^a,Xa) = " [{C + DM)-'D]^^ fl - [{C + D D] fahe''^ 

- [iC + DM)-'D].^^ faf-ae^'-^ - [{C + DM)-'D]^^ f? + ^Jah (10) 
Integration over the Graph may be decomposed as integration over the bonds : L.^.^ ^ dx Q{x, x) = 



X^(a) /o" "i^a Qa,aixa, Xa)- The integration of product of functions fa and fa is given by eqs. (71 
hence 



72) 



G(7) = / dxQ{x,x) = Ti{{C + DMy^Dd^M} -d. 

J Graph 



(a) 



H-faiL)] 



(11) 



where the trace runs over the 2B arc indices. Up to now, nothing has be assumed on the 
nature of the boundary conditions. An expansion of the trace provides a trace formula for 
/craph Q{x, x), where the expansion is interpreted as sum of contributions over cycles (orbits) 
in the graph. 

We now suppose that the matrices C and D do not depend on the spectral parameter 7. Per- 
forming and integration In 5(7) = f"' d^y' G{'~f'), we obtain the first important result of Desbois : 




(12) 



1. In one dimension, a dependence of the vector potential in the coordinate may always be removed by a 
convenient gauge transformation. 



4 



where the product runs over the B bonds. This provides a general expression of the spectral 
determinant in terms of 2B x 2i?-matrices coupling arcs. We repeat that the matrix M couples 
reversed arcs and contains local information related to the potential on each bound. On the 
other hand, the matrices C and D, characterising which arcs arrive and issue from the vertices, 
encode the information on the topology of the graph. Note that the structure det(C + DM) 
was also obtained in [19] in the absence of the potential where the determinant det'(— A) was 
analysed (the prime indicates exclusion of zero mode). When V{x) = the matrix M takes the 



form Ma^b = -y/l{5a,b COth. - 5, 



a,b sinhy^ia 



The prefactor in Eq. (12) has been chosen arbitrarily and for convenience for the following. 



We will come back to the question of the prefactor in the section |6j 

Limit 7 — )• CX3 For large positive spectral parameter we can ignore the potential, therefore 
/^(O) ~ s-'^d fai^a) — — ^/^e"^'". Hence the matrix M becomes diagonal : Ma^b ~ 

—Safis/l we obtain 

5(7) ~ i-lf 2-^ -f-^/^e^^ det{C -J^D) (13) 

7— >+oo 

where ^ = la is the "total length" of the graph. 

Example Let us consider a ring of perimeter L pierced by a magnetic flux 9 (a graph in which 
the number of vertices can be reduced to ^ = 1 with B = 1 bond) . We denote by 1 and 1 the 
two reversed arcs. Let us consider boundary conditions described by the matrices given below 



by Eq. (29). We obtain 



A-/((0)-/i(0) + 2/((L) cosg 
-f[{L) 



(14) 



It is interesting to point that this formula has found a practical physical application for a 
potential of the form V{x) cx x{L — x) in the context of the study of decoherence by electron- 
electron interaction in a phase coherent metallic ring [l3] (hence /i(x) is a Hermite function and 
fi{x) = fi{L — x) due to the property V{x) = V{L — x)). 
The case V{x) = is analysed in section [6] 

3.2 Arc determinant (2) 

The above derivation has privileged the use of variables Qa giving the value of the Green 
function at the vertices. Another natural choice is to deal with (covariant) derivatives of the 
Green function at the vertices. We write 



5a(max (xa, Va)) ^.^(max {Xa, Va)) 



(15) 



where the functions Qa and Qa are solutions of the Schrodinger equation on the bond (a) for fixed 
values of the derivative at the boundaries, 5^(0) — 1 and g'a{la) = (see appendix [A|) . Following 
the lines of the previous subsection, we obtain 

/ : \ 



Q{o) = Ng\o) 





gb{yb)e-'^>'y^ 



(16) 
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where we have introduced the arc matrix 

Na,b = 5a,bgai0) + <7a(/a) 6^'-^ . (17) 

Imposing the bomidary conditions at the vertices we obtain the values Reinjecting these 



expressions in Eq. (15), integration over the coordinate and summation over bonds give another 



general trace formula for the Green function 



G(7)= [ dxg{x,x) = TT{{CN + Dr^Cd^N} -d^Y.ln[-ga{la)]. 

JGraph 

For 7-independent matrices C and D, integration over the spectral parameter leads to 

9a\''a) 



(18) 



(a) 



(19) 



where the 7-independent prefactor was chosen in order to match with (12). This expression can 



be directly related to (12) by using the relations (demonstrated in the appendix [A[) 

-1 1 _Tjf'aila 



N = M 



and det M 



1 

detiV 



n 

(a) 



9a{la) 



(20) 



Example We consider again the very simple case of a ring, but this time we choose to consider 
boundary conditions of the type described below by Eq. (|40|). We obtain 



/^-gi(0)-gi(0)-2gi(L) cosg 
-9i{L) 



In the absence of potential we obtain 

5(7) = /x-y/7sinh ^/^L + 2 (cosh ^/^L + cos ( 



The limit — )• 00 gives the Neumann determinant of the bond : 5(7) — 

/^^7n:r=i(i + (#)^)- 

3.3 Arc determinant (3), scattering matrices and (^-function 



(21) 
(22) 

/x^sinh^L = 



The result (12) can be reorganised more conveniently by introducing the arc matrix 

i2 = (^ + M)(V7-M)~^ (23) 

i.e. M = ^{R + 12b)-^{R - I2B) = ^/7(« - ^2b){R + l2i?)"\ where I2B is the 2B x 2B 
identity matrix. We denote its matrix elements 



Ra,b = Sa,b ra+S^ltae 



Ida 



(24) 



These matrix elements have a clear physical meaning : and ta are analytic continuations (to 
negative energies E — )• —7) of reflection and transmission probability amplitudes through the 
potential Va{x) (see appendix [A|) : R is therefore the analytic conitnuation of the bond scattering 
matrix. We also introduce the analytic continuation of the vertex scattering matrix (scattering 
matrix interpretation is discussed in Refs. [Ml SB) 



Q 



C + ^D)-\C + ^D) 



(25) 
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(for instance we can check that Q is unitary for = —ik, reflecting current conservation at 
the vertex) . We obtain another interesting result of Desbois [10] : 



||/a('a) det{l2B + R) 

(a) 



(26) 



Let us describe the structure of this result : the product over bonds Yl(^a) f'ai^a) deter- 
minant det(l2_B + R) involve (local) information about the potential on the bonds (recall that 



R couples an arc to itself and its reversed arc, only). This is even more clear from Eq. (78) that 
leads to 

n[-/a(^'^)] det(l2B + R) = 2^^^/'llt. = (47)^/2 [n Raraf" , (27) 

(a) (a) 

where the first product runs over the B bonds and the last one over the 2B arcs. Now let us 
consider the determinant det(C— ^ D) : organising the arcs by gathering arcs issuing from the 
same vertex, the matrices C and D take some block diagonal structure, hence det(C — ^/j D) 
encodes some (local) informations about the vertices. The last determinant det(l2B — QR) 
mixes informations about potential and the connection of arcs to vertices in a nontrivial way, 
and characterises the (global) information about the topology. 

This last part of the spectral determinant vanishes on the spectrum of — + V{x), i.e. for 
7 = —En- Because Q and R have the meaning of scattering matrices, the equation det(l2s — 
QR) = may be understood as a quantisation condition a la Bohr-Sommerfeld. 

It is well known that the determinant det(l2_B — QR) can be expanded in terms of primi- 
tive cycle contributions, with the structure of a ("-function. We recall that a cycle (a periodic 
orbit) is the equivalence class of all ordered sets of arcs C = (ai, 02, • • • , an) identical by cyclic 
permutations and such that Vi € {1, • • • ,n} end(ai_i)=begining(aj) (with oq = an)- An orbit is 
said primitive, and denoted C, if it cannot be decomposed as a repetition of a smaller orbit. We 
define the weight of the orbit C as v{C) b(C) = {QR)aia2iQR)a2a3 ■ ■ ■ {QR) 

anai where we have 

identified the parts related to the scattering by the vertices v{C) = QaibiQa2b2 ' " Qa„bn and 
the scattering by the bonds b{C) = 02-^^2(13 ' ' ' Rbnan respectively (note that G {a^, ai}). 
The determinant may be rewritten as an infinite product over primitive orbits as 

det(l2ij -QR) = 1[{1- v{C) 6(C)) (28) 
c 

(see for example Refs. \37\ [T] and references therein). This relation emphasizes that the spectral 
determinant may be interpreted as a (-function (or L- function) for primitive orbits of graphs, 
(-functions are powerful tools of particular importance, in number theory and graph theory 
for example, since they play the role of generating functions for primitive elements. The most 
famous (-function is the Riemann (-function ((s)~^ = riprimepC-'- ~ P~^) giving access to the 
distribution of prime numbers. 



To close this section, let us emphasize that the representations ( 12|19|26 ) are remarkable in 



the sense that, despite the Laplace operator acts on a space of infinite dimension, its spectral 
determinant may be related to determinants of finite size matrices : all above expressions have 
involved arc matrices of size 2B x 2B- We will now see that we can express the spectral 
determinant in terms of ^ x y vertex matrices in general smaller]^ 

2. The spectral determinant can only be expressed in term of vertex matrices when boundary conditions at 
the vertices are such that it is possible to assign a unique variable to each vertex. This can only be done for 
"permutation invariant" boundary conditions ilQ, . Two such particular cases are discussed in the two following 
sections. 
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4 Continuous boundary conditions 



4.1 Roth's trace formula 



Let us consider the important case of the Laplace operator acting on functions that are con- 
tinuous at the vertices : V'a/3(0) = ''Pa for aU vertices f3 neighbours of a, and aai3(Dx'ip)ai3{0) = 
Xaipa, where the connectivity matrix ensures that the sum runs over vertices neighbours of a. 
The parameter must be chosen real in order to ensure self adjointness of the Schrodinger 
operator ; it can be understood as the weight of a (5 potential at the vertex (these boundary 
conditions are sometimes denoted as "(5-coupling" [12^ [T3] ) . In the appropriate basis where arcs 
are gathered by vertices from which they issue, the matrices C and D have block diagonal struc- 
tures (with V blocks a = 1, ■ ■ ■ V) where the x nia block a correspond to the vertex. A set 
of 7-independent matrices C and D corresponding to (5-couplings are made of x rria blocks 



-Xa 

-1 1 
-1 1 



V -1 



\ 




1 / 



and 



/I 1 





V 



1 1 \ 



... 

■•• 0/ 



(29) 



It will be useful to remark that det(C(") - ^D^")) = -(A Q, + niQ.^^). The vertex "scattering 
matrix" for the vertex g(") = (-C(") + ^D("))-1(C(") + takes the form [I] 



Q 



(a) 



/I 
1 



V 1 



1\ 

1 



V 



<-m.a ) 



(30) 



— for visiting the vertex a and 
oo corresponds to Dirichlet boundary 



(for 7 = we recognise the weights introduced by Roth 

— 1 for a reflection on it). Note that the limit A^ = 

conditions ipi{0) = 0. 

Moreover let us consider the case V{x) = 0, then = and ta = e~v^'" (note that this latter 
is the analytic continuation of the transmission probability amplitude through a free interval 
ta = e'*^'" for = \J —k"^ 



iO+). Starting from (26) and using (78) we get 



5(7) 



e^^Jj(l-a(C)e- 



V7^(C)+i0(C) 



(31) 



det(l2s-Q/?) 



where the first product runs over the V vertices. L = '}2i{a) is the total length of the graph. 
The weights of the orbit C = (ai, a2, • • • , a„) is now denoted v{C) —?- a{C) = QaiaiQaia-^ ■ ■ ■ Qa„ai 



for Q given by gathering the blocks (30). The notations £{C) and 9{C) designate the length of 



the orbit and the magnetic flux enclosed by it, respectively. Thanks to the simple structure 
of the matrix R for V{x) = 0, the bond scattering part of the weight of the orbit is simply 

5(c) Q-vim+im, 

If we choose moreover vertex scattering with Aq = 0, the weights a(C) become 7-independent, 
what allows a simple Laplace inversion of the derivative ^ In ^(7) = dt Z{t) e"''*. We recover 
the Roth's trace formula 



^ y-B 1 

2^/7rt 2 2Vvrt ^ 



£(C)a(C)e-T^ 



hie{C) 



for Ar 



OVa 



(32) 
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where the sum now runs over all periodic orbits C, where C is the primitive orbit related to the 
orbit C (see Refs. [Mill] for more details). 

4.2 Vertex determinant 

Another interesting expression of the spectral determinant in terms of x l^-matrix coupling 
vertices may be obtained by a construction anologous to the one given above. The asumption 
that the functions are continuous at the vertices allows us to deal with vertex variables (see 
figure [l}b) : 



QiAai3{Xai3-yc,l3) 



+ 6, 



w, 



/a/3 (max {xap^Vap)) //3a (max {xp 

)) (33) 



where Ga is the value of the Green function at the vertex. The boundary condition at vertex ^ 
takes the form 



Mf,u Gu = S^a fapiyap) e"^^"''^"^ + fpa{ypa) e"'^""?'''" 

V 

where the V matrices coupling vertices is given by : 

A^a/3 = <5a/3 - ^ aa/./aM(0)^ + f'l3ai^»f}) e"'^"" ■ 



(34) 



(35) 



We can now obtain the value of the Green function at the vertices and reinject these expressions 



in (33). Integration of 



^a/3,a/3(a;a/3,a;a/3) = {M ^)^^ fal3{Xal3f + W^"^ {M ^)„^ + e '^"^ {M 1 f aliiXa^) f Pa{x ^a) 



+ (-^ ^)/3/3 fpai^liaf + ^^^Jaf}{Xap) fpa{xpa) 



w, 



a(3 



can be performed along the same lines as before. One obtains [H] 

S{l) = ll^-f'ap{lap)]-' detM 



(36) 



(37) 



where the product runs over the bonds of the graph. Eq. (37) is less general than Eq. (12) since 



it corresponds to a particular choice of boundary conditions at the vertices, however it involves 
the information about the graph in a more compact manner in the sense that the V xV matrix 
M is smaller than the 2B x 2B matrix C + DM. 

In the absence of a potential, V{x) = 0, we recover Pascaud & Montambaux' result 



{aP) 



^/7 



for 



M 



al3 



Sal3 [ ^a + y/l'Y Qafi COth 



a^ I fla/3 



sinh ^lai3 



(38) 



(39) 



Note the alternative derivation of (38) and (37) with the path integral 11] and [9]. It is worth 



pointing that the full spectrum of the Laplace operator is not always given by det = : it 
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may occur in some particular cases that eigenstates of the Schrodinger operator vanishes at all 
vertices ipa = for ip{x) ^ 0. Some examples are discussed in details in Ref. [IQ] (this is related 
to the phenomenon known as "bound state in the continuum" in the scattering situation, for 
graphs with some infinitly long wires \A5\ l38] ) . 

Finally, note that the 7-independent prefactor in (37) is chosen in order to match with 
formulae (12) and (26) : one can easily check that the 7 — )• 00 behaviours coincide. Note that, 
in Ref. [Ij, a direct relation between (38) and (31 ) was established, without a posteriori matching 
procedure. 



5 Continuous derivative boundary conditions 



5.1 Trace formula 

Another interesting case worth to be discussed is the case of Schrodinger operator acting 
on functions whose derivative are continuous at the vertices {T)xil^)ai3{0) = V'q fo^ vertices 
/3 neighbours of a, and aa^V'a/3(0) = fJ-a^Pa- This choice is denoted as "(5^-coupling" in 
Refs. [T2|[T3]. A set of 7-independent boundary matrices at vertex a are in this case 



( 1 1 





V 



1 1 \ 



... 

•• 0/ 



and 




110 
1 1 



\ -1 



\ 



1 / 



(40) 



that gives det(C(") - ^Z^^")) = (-l)"'"+S'""'^^(Aia + "^a/Vr)- F^om ([25j), we get the vertex 
"scattering matrix" for the vertex a : 



Q 



(a) 



(\ 1 ... 1\ 
1 1 ... 1 



m„ + /Uai/7 



(41) 



\\ 1 ... 1/ 

therefore the limit fia = 00 corresponds to Neumann boundary conditions il^i{0) = 0. 



Starting again from Eq. (26), a little bit of algebra gives the spectral determinant in the 
absence of potential 



-v/7^{C)+i0(C)j 



(42) 



where the weights v{C) — t- /3(C) are computed with matrix elements of (41). It is interesting 
to point that these weights may be simply related to the weights a(C) of the previous section 
thanks to the substitution Aq, — )• 7/^0 and adding a sign (|_i)#arcsof c_ 

If we moreover set /i^ = we may obtain the trace formula for the partition function 



C B-V 
2v/7ft 2 2Vvrt 



^^^(C)/3(C)e- 



for /liq, = V a 



(43) 



that generalises the result of Roth to the case of the Laplace operator acting on functions whose 
derivative is continuous at the vertices. Compare to the case of functions continuous at the 



vertices discussed in the previous section, Eq. (32), only the constant term and the weights of 
the orbits have changed in sign. 



10 



5.2 Vertex determinant 



The expression of the spectral determinant may also be obtained by using vertex vari- 
ables 



(a/3) 



with 



(44) 



(45) 



where the functions Qa were introduced above (and in the appendix [A|) . In the absence of the 
potential, V{x) = 0, we obtain 

(46) 



Sil) = Yi \/7sinh^/„/3 detTV 
where A/"a/3 = Sa/s [fJ-a + ^ "-aiy coth ^lau) + Oq/? 



^sinh y^l 

a.p 



6 From the trace formula to the (^-regularised determinant 

The question of the prefactor of the spectral determinant is the subject of the present section. 
It is important to fix the 7-independent prefactor : in all expressions given in the previous 



sections, (12), (19), (26), (31), (38), (42) and (44), the prefactor has been chosen arbitrarily. 



nevertheless in a way that all expressions match when varying continuously boundary conditions. 
However this choice is a priori not related to a particular regularisation and has only been made 
for convenience. 



6.1 Functions continuous at the vertices 

Let us first consider the simplest case of the Laplace operator acting on functions continuous 
at the vertices in the absence of potential. Moreover we set the parameters Aq = 0, when the 
Roth's trace formula for the partition function holds. Mellin transform of the Roth's partition 
function (32) reads 

Xs-1/2) . V-B 

7 



r(-i/2)r(s) ' 2 
+ -r'^"' E 'iCMcW"^' (^) '<i-. (vim) (47) 

where Kiy(z) is the MacDonald function (modified Bessel function) [iSj. First and third terms are 
of the form p^/i(-s) for a function h{s) regular for s — )• 0, therefore we can use [f^/i('S)] q = h{0). 

The analysis of the third term requires the formula [^^]^=i/2 = -v^e^Ei(-2z) Using 



and Ki/2{z) = yji ^ ^ we get 

^(0,T) = -v/7£ - ^ nl^^e-V^'C^'+'^W (48) 

Sum over orbits can be decomposed as a sum over primitive orbits and their repetitions = 
J2cT^n=v Using i{C) = ne{C), 9{C) = n9{C) and a{C) = a(C)". Therefore 

dC,„ , V-B 



ds 



(0,7) = -^C- "—^ ln7 - ^In 1 - a(C)e-v^^('^)+'^('^) (49) 
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that is 



7 



^ev^^J|(l-a(C)e 



-^£{c)+ie{c) 



7 2 



ev^^ det(l2B - QR) 



(50) 



We recognise part of the expression (31) (see also Eq. (65) of Ref. [Ij). We have related above 
3r representation in 1 

2 sinh y^/a^ det Ai 



Eq. (31) to another representation in terms of a vertex matrix, Eq. (38), hence 



5^(7) = n 



{a/3) 



\/7 Ua^c iUa^a) 



S{j) for Vix) = 0k Xc 



Va 



(51) 



The previous formula may be easily generalised to the case of the Schrodinger equation 
for continuous boundary condition with arbitrary parameters Aq,. Using the observation that 
various regularisations only differ by a 7-independent prefactor, we conclude that 5^(7) = 
2^ / (Yla f^a) S{'y) also holds in the presence of the potential 



det 



(52) 



This is precisely Friedlander's result |17j (c/. Ref. [44j, where the correspondence of notations 
is discussed). The product 11(0/3) ["^//q/jI^o/?)] the Dirichlet determinant of the graph (de- 
terminant when Dirichlet boundary conditions are chosen at all vertices that make all wires 
independent). 

Example 1 : wire Let us consider a wire of length L {B = 1 bond and V = 2 vertices) with 
no potential V{x) = 0. Matrices C and D for boundary conditions of type ([29j) are 



C 



-Ai 






-A2 



and D 



1 
1 



where Ai and A2 describe the boundary conditions at the two vertices. Eq. (61) gives 



^wire(7) = 2^/7 sinh + (Ai + A2)2 cosh + A1A2 



2 sinh ^/^L 
V7 



(53) 



(54) 



This result can be obtained straightforwardly from (|52|). 
For Ai 



A, 



3(7) = 2/7 sinh 



we recover the Neumann determinant S^jj-ev 

(t)^ with n G N are the eigenvalues of the Laplace oper- 



27inr=i(l + ii;) where : 
ator on (0, L) for Neumann boundary conditions. 

Note that the second term 2coshy^L corresponds to the Neumann/Dirichlet determinant 
(Neumann at x = and Dirichlet at x = L or vice versa) retained for Ai = and A2 = 00 (with 
A1A2 = 0). 

The last term ^ ^'"^^^^^ is the Dirichlet determinant (Ai = A2 = 00). 

Example 2 : ring We consider a ring of perimeter L (B = 1 bond and V = 1 vertex) with 
no potential V{x) = and pierced by a flux 9. Matrices C and D for boundary conditions of 



type (29) are 



C 



-A 
-1 1 



and D 



1 1 




(55) 
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therefore (61) gives 



5ring(7) = 2 ( cosh - COS + A 



sinh 

V7 



(56) 



It may also be obtained from (52) knowing how a loop contributes to M. (see Ref. p[]). 

[l + where 



For A = we obtain 5'^ing(7) = 2(cosh^/7L — cos 6*) 
^2n|^^2 ^-^j^ n G Z is the spectrum of the ring. 



Example 3 : star graph We consider a star graph with B arms (then V = B ^X) with no 
potential V{x) = 0. We introduce a parameter A at the central vertex and set all other boundary 
parameters equal to zero, Aq = 0. In this case it is more easy to use (51). The B wires are 



"dead arms" and may be taken into account very simply through the following rule [32] : a dead 
arm issuing from the vertex a gives a contribution ^77 tanh y^la that should be added to the 
matrix element A4aa corresponding to the graph from which the dead arm is removed ; in the 



product over bonds of Eq. (|51|) we replace by cosh ^la- This simple rule can be easily 



^ + \/7Ef=itanh^/a 



recovered using the path integral formalism introduced in Ref. [T]. It allows to treat the star 
graph as a graph with one vertex and B dead arms. Therefore M 
and the spectral determinant is given by 



'S'ftar (7) 



2 

^B 



B ^ / 
- Jl cosh A + 



a=l 



B 

7 tanh 

a=l 



(57) 



If we consider the case A = and take the limit 7 — t- we recover the result of Ref. [19] 

SiM -^7 = 7 det'(-A) 
where det'(— A) is the determinant when the zero mode contribution has been excluded. 



(58) 



6.2 Functions with derivative continuous at the vertices 

The reader has remarked that there is very little difference between the Roth trace formula 



(32) and its generalisation (43) : the second term proportional to F — i? is changed in sign 
and the definition of the weights are changed : when Aq, = and /Uq = we have simply 
/3(C) = (— l)*'^'^'^'^°^''a(C). We can therefore perform a similar calculation and get in the case 
V{x) = and Ha = : 



5^(7) = 7^e^^ n (1 - /5(C) e-v^^('^)+'^('^)) (59) 
c 

By a similar continuity argument as in the previous subsection we get for V{x) 7^ and ^q, 7^ : 

(60) 



5«(7) = n 



detM 



{al3) 



the first term (product over the bonds) coincides now with the Neumann determinant of the 
graph (product of spectral deterimants for all disconnected wires with Neumann boundary con- 
ditions). This latter expression hence provides some generalisation of Friedlander's result to 
another kind of boundary conditions. 
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7 Conclusion 



Using the Roth's trace formula for the partition function and its generalisation, we have 
obtained the expression of the (^-regularised spectral determinant of the Schrodinger operator 



acting on functions continuous at the vertices, Eq. (52), or functions whose derivative is contin 



uous, Eq. (60). It is worth emphasizing that in these formulae, the 7- independent prefactor is 
now fully determined by some properties of the Graph : i.e. the ("-regularisation fixes a priori 



the prefactor of the spectral determinant, whereas in formulae (12), (19), (26), (31), (38), (42) 



and (44) it was chosen arbitrarily a posteriori. 

We have summarised the relations between the various expressions of the spectral determi- 
nant given in the article in figure [2] 

The work of Desbois [10\ on the case of general boundary conditions, in particular equations 
( T2|l9|26 ) , allows in principle to go continuously from one to the other of these two situations by 
modifying continuously the boundary condition matrices C and D. Therefore, by continuity this 
leads us to conjecture that the (^-regularised spectral determinant of the Schrodinger operator 
for the most general boundary conditions is given by 



(-1) 

n 



V 



n 



det{C + DM) 



det(7-i/^C7 - j'^/^D) 



det(l2B - QR) 



(61) 



where runs over the V vertices, 0(0) '-'^^'^ ^ bonds and Yl^ over the 2B arcs. It would 
be interesting to provide a direct proof of this formula. 



Eq. |[T2j i 

S ocdct{C + DM)\ 

X"' 



Eq. | |37| 

S ocdetM 



\Eq. ^ \ 
"\ S ocdct{CN + D)\ 



Eq. |[26]l 

S oc det(l - QR) 



Jpix) continuous 
at vertices 



1/ , 



Eq. 



Trace formula 
for Z{t) (Roth) 
Eq. JMI 



5<; 

Eq. |52l 



X 



^'{x) continuous 
at vertices 



Eq. l|44| 

S oc detAf 



Eq. |p2l 



I 





Trace formula 




for Z(t) 




Eq. (pS) 




1 










Eq. (|60| 






/ 



Conjecture : Eq. ||61|i 



Figure 2: Relations between the several expressions derived in the article. The dotted green boxes 
indicate results valid in the most general case whereas blue boxes are valid in particular cases. 
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A Bond scattering 

Let us consider the Schrodinger equation on the bond (a) for an energy —7 = E = : 



d' 



2 



^ - d^! + 



V^(Xa)=0 iorXaG{0,la) (62) 
We introduce three useful basis of solutions of this differential equation. 



Basis 1 We denote by (j)a{xa) the scattering state with an incoming wave from the left, with 
(l^aixa < 0) = e^'^^" + raQ~^^^°- and <t)a{xa > la) = t^e^'^^^""'") = tae~^'^^^, where r^, ta are 
reflection and transmission amplitudes through the potential : 

0a(O) = l + ra and cpaiL) = ta . (63) 

The scattering state incoming from the right is naturally denoted by (pa{xa)- 

A well known sum rule relating integral of the square modulus of the wave function (i.e. 
local density of states in the scattering region) to scattering matrix is the Krein-Friedel sum rule 
|16 1 127 ^ 136]. Generalisations of this sum rule in the context of metric graphs have been obtained in 
Refs. jiniSI]. For E = k"^ = —7 > 0, let us change notation and introduce '^'^^\x) = ~^j^^4>a{x) 

and ^^^\x) = ^^j=^4>a{la — x), the left and right stationary scattering states, respectively. We 
can deduce from the local version of the Krein-Friedel sum rule demonstrated in Ref. [41] for 
graphs that 



dx [^("\x)Y^^f'Hx) = — [y^^ + j ^ with a, /3 G {L, R} (64) 



where =5^ is the 2x2 scattering matrix describing the scattering by the potential on the bond 



We obtain 



dx \Uxr = -2iVE + + (66) 



J^^dxMxrMla-x) = -2iVEl^r:^ + t:^ + '^^ . (67) 

Basis 2 Let us introduce the solution fa{xa) satisfying the boundary conditions 

/a(0) = l and /a(/a) = 0. (68) 

It follows that fa{x) is a real function for 7 G M. For example, when V{x) = 0, the function 
reads fa{x) = '^'"'iiA^f • Another independent solution of this differential equation is naturally 
denoted fa{la — Xa) = fa{xa) and take the values for Xa = and 1 for Xa = la (do not confuse 
the two functions fa and fa with the components of a scalar function). The Wronskian of the 
two solutions is 

Wa = W[fa{x), f-a{la " x)] = -f'M = - f'Sa) • (69) 
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Sum rules can be obtained as follows [S] : we remark that d.yfa{x) is solution of the differential 
equation [j- + Va{x)]djfa{x) = -fa{x) for the boundary conditions djfa{0) = d^fa{la) = 0. 
Integration straightforwardly gives 



d^ifa{x) 
We deduce 



1 



fa{x) / dx' fa{x')f-a{la ~ x') + Uk " x) / dx' /i(x') 



dxfa{x)faila - x) = d^faiQ ■ 



(70) 
(71) 

(72) 



They replace the Krein-Priedel sum rule like formulae given above for scattering states. 

One can establish the relation between the two basis of solutions {(paixa), (j^aixa)} and 

{fa{Xa), f-a{Xa)} [39j : 

X {'^ + ra)(l)a{^a) - ta(l)a{Xa) .„„x 



{l+ra){l + ra)-tata 



The two following relations follow 



f (n\ = - r (i-^a)(i + ^a) + tata 



faiQ 



2ta 

{l + ra){l + ra) -tata 



(74) 

(75) 



where we have performed some analytic continuation to negative energies 7 = -k^ - iO+. Let 
us introduce the 2x2 block in the matrix M related to the two arcs a and a : 



M, 



(a) 



(76) 



The relations ( 74|75 ) with (23) immediately shows that the 2x2 block in the matrix R related 
to the two arcs a and a indeed encodes the reflection and transmission amplitudes : 



(a) 



ta e 



tae 



that is Ra,b = Sa,b ra + S^^tae 



(77) 



Hence demonstrating that R is the "bond scattering matrix" (or more precisely its analytic 
conitnuation to energy E — )• —7). We obtain the useful relation 



f'M det(l2 + %)) = 2V7ta 



(78) 



(we recall that ta = ta follows from the one-dimensional character of the wire : it can be 
demonstrated by computing the Wronski determinant of the two scattering states at the two 
edges of the interval) . 



Basis 3 Another useful set of solutions are the functions ga{xa) and ga{xa) whose derivative 
take values 1 and at the two ends of the interval : 



9'aiO) = l and g'aila) = 0. 



(79) 



For example, when V{x) = 0, the function reads ga{x) - 
nant is given by 

= W[ga{x),ga{la " x)] = -ga{la) = -gM • 



^°J^£!h^U - ^^"^ Wronski determi- 



(80) 
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Following the method of Ref. [H] and recalled above one can construct the derivative of the 
function with respect to the spectral parameter d^ga{x) = [Oaix) Jq dx' ga{x')ga{la — x') + 

g-a{la - x) Jl^ dx' glix')] . Therefore 



Ix 

'•la 







dx gaixY = d^gaiO) (81) 



la 

dxga{x)ga{la - x) = d^gaiD (82) 
It is also interesting to relate this basis of solutions to the two other basis : 

g ^ ^ 1 {1 - ra)(j)a{Xa) + ta(j)a{Xa) ^^^^ 

^ (1 - Ta)(l - ra) - tata 



from which we obtain 



1 (1 + ?'a)(l - r-a) + tat-a , . 

5a0=-^7^ 84 

^ (1 - ra)(l - r-a) - tata 



We also easily demonstrate that 



from which 



fMfa{Xa)+fa{la)h{x-a) 



^"^^^ mfm-ma)ma) ^^^^ 

'^^^"^ ^ imrnf^imim ^^^^ 



The two matrices introduced above are hence simply related by 

Ida 



AT =( aM.^ ^^(''^^ ^ " ) = M-( . (89) 



Note also that detM(a) = faiD/gaiD- 



B (^-functions for combinatorial graphs 

The study of ^-functions (also denoted L-functions in the mathematical literature) and trace 
formulae in combinatorial graphs has attracted a lot of interest [22l EOj EH El |37l [2] (they found 
recently some application in the context of quantum chaos [29t I30| ) ; a brief recent review on 
mathematical aspects may be found in the introduction of Ref. |28j . We show in this appendix 
that the general trace formula for combinatorial graphs (the Bartholdi's formula [2] generalising 
Bass [3] and Ihara [22j trace formulae) may be deduced from the trace formula obtained in 
Ref. UJ for metric graphs. This latter is itself a particular case of the general trace formula 
obtained by Desbois [TU] recalled in section [sj 

We consider the trace formula for the Laplace operator acting on functions continuous at 



the vertices. Our starting point is the equality, Eqs. (26) and (37) 



^_^^^detg7_^7l^ det(l2B - QR) = detM (90) 
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We recall that in this case (—1)^ det(C — ^/^ D) = W^{\a + ^a^ff) 



Note that (90) is related to the spectral determinant of the metric graph when bomidary 



condition matrices C and D are 7 independent. However, having established the relation between 



the vertex determinant and the arc determinant, Eq. (90) also holds when the parameters 
depend on 7 (the case considered below). In this case, Eq. (90) is however not anymore related 
to the spectral determinant of the metric graph. 

Let us consider the case with no potential V{x) = and set all the lengths of the wires 
equal : la = l'i a. For convenience, we introduce the notation 



wu, (91) 



where u and w are two real or complex numbers. The bond scattering matrix elements are 
Rij = wu6^ j and therefore det(l2B + R) = [I — (uw)"^]^ . We consider the case where the 
boundary condition parameters Aq, are related to the valencies of the vertex by the relation 

ma + — = 2w. (92) 

The matrix QR has elements : {QR)ij = n if the arc j terminates at the vertex from which 
issues the arc i (with j ^ i) ; {QR)^^ = u{\ — w) ; all other matrix elements are zero. We write 
(notations reminiscent of those of Refs. [21 128| ) 

QR = u{^ -w^) (93) 

where = 1 if end(j) =begining(i), = otherwise ; the matrix ^ couples reversed arcs 



On the other hand, from (39) we find 



Map = 7 — ^ Uai3 [w + {rria - w){uwf] - uwaap] (94) 

Using that (-l)^det(C - ^ D) = {2w^Y 

we deduce the Bartoldi's formula for the zeta 

function [21 [JOl |28l |29] 

Z{u, w)-^ = n " " ^)"''^^^ ^^^^^) (95) 
c 

= det(l2B - - w J)) = (1 - {uvof)^-"^ det(li/(l - {uwf) -uA^ wu^y) , (96) 

where i{C) is the length of the orbit (number of wires visited) and n^{C) the number of reflections 
of the orbit on vertices. A denotes the adjacency matrix with matrix elements a^^. The matrix 
y is the diagonal matrix gathering the valencies : 3^q,^ = maSai3- This formula was used in 
Ref. [\.0\ [7] as a generating function for counting of the orbits with finite number of reflections. 
It has also found some interesting application in the context of quantum chaos in combinatorial 
graphs ^29^ i30j . 

Note that combinatorial graphs are related to metric graphs when all lengths are equal 
and for permutation invariant boundary conditions at the vertices. It follows that only two 
weights can be attributed to an orbit passing at a vertex (u for no reflection and u(l — w) for a 



reflection) and therefore (97) is the most general trace formula for combinatorial graphs (up to 
the straightforward addition of magnetic fluxes in matrices A and ^) . As a consequence all the 
trace formulae obtained for metric graphs with permutation invariant boundary conditions at 
the vertices (we have consider the particular case of continous boundary conditions here) would 
lead to the Bartholdi formula when lengths are taken to be equal and potential vanishes. 
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For w = 1 we recover the Bass formula O [37] (this was also discussed in Refs. [TUl E]) 



for the ("-function over the backtrack- less orbits Cb (orbits with no reflection). The Bass formula 
generalises the Ihara formula [22] demonstrated for regular graphs. 
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